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Abstract 



Controversy concerning the pairing symmetry of high-T c materials has 
motivated an interest in those measureable properties of superconductors for 
which qualitative differences exist between the s-wave and d-wave cases. We 
report on a comparison between the microscopic electronic properties of d- 
wave and s-wave superconductors in the mixed state. Our study is based on 
self-consistent numerical solutions of the mean-field Bogoliubov-de Gennes 
equations for phenomenological BCS models which have s-wave and d-wave 
condensates in the absence of a magnetic field. We discuss differences be- 
tween the s-wave and the d-wave local density-of-states, both near and away 
from vortex cores. Experimental implications for both scanning-tunneling- 
microscopy measurements and specific heat measurements are discussed. 

PACS: 74.60.Ec, 74.72.-h 
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Since shortly after the discovery of high-temperature superconductors (HTSC), there has 
been great interest in determining the pairing symmetry of the order parameter |IJ . In the 
absence of disorder, low temperature electronic properties in the Meissner state of a d-wave 
superconductor differ qualitatively from those of a conventional s-wave superconconductor 
because of the existence of nodes in the gap function. These differences can in principle 
be used to identify the pairing symmetry, although strong anisotropy and the complicated 
nature of the materials have conspired to make conclusive experiments difficult. (Recent 
work is strongly suggestive of d x 2_ y 2 pairing.) It is also of interest to study differences be- 
tween the mixed-states of d-wave and s-wave type II superconductors. In the mixed-state 
magnetic flux will penetrate the superconductor and form an Abrikosov vortex lattice. Low- 
lying quasiparticle excitations will then exist for both pairings, although in the conventional 
case they must be bound to the vortex core where the order parameter vanishes. The exis- 
tence of bound quasiparticle states was first predicted by Caroli, de Gennes and Matricon 
when they studied an isolated vortex line in a conventional superconductor. Experi- 
mentally, these quasiparticles have been observed in scanning-tunneling-microscopy (STM) 
measurements ||. For the d-wave case, progress has recently been made on both experimen- 
tal and theoretical fronts. Volovik has used semiclassical approximations to calculate 
the density-of-states (DOS) at the Fermi energy N(0) for the mixed state of a d x 2_ y 2 su- 
perconductor in a weak magnetic field H -C H c2 . He found a finite N(0) in the absence 
of disorder proportional to H 1 ^ 2 compared to the H 1 behavior expected for conventional 
superconductors in the same approximation. This prediction appears to be in accord with 
recent measurements of the magnetic field dependence of the low-temperature specific heat 
H|7| in high T c materials. However, the short coherence length of high T c materials raises 
some uncertainty about the detailed applicability of a semi-classical analysis and motivates 
a fully microscopic study of the same problem. In this Rapid Communication we report on 
such a study. 

Application of microscopic mean-field-theory to inhomogeneous states of superconduc- 
tors gives rise to the Bogoliubov-de Gennes (BdG) equations |§. Motivated by the STM 
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experiments of Hess 0, numerical solutions of the BdG equations have been obtained for 
both continuum ||[10| and lattice [TTJ models of a superconductor containing an isolated 



vortex. This work has recently been generalized to the case of isolated vortex in a d-wave 
superconductor |12|. According to Volovik, the DOS in the mixed state of a d-wave su- 
perconductor is dependent on the typical distance between vortices so that for the present 
study it is necessary to solve the BdG equations for the vortex-lattice state of a d-wave 
superconductor. 

To model decoupled C11O2 layers we consider single-band Hamiltonians on a 2D square 
lattice with nearest neighbour hopping and both on-site and nearest-neighbor interactions: 

H = H° + H', (la) 

h° = - (tij ]^ + ^ c *W) -E^> ( lb ) 

<ij><7 icr 

H' = U ^ n^hn + — h ia h ja >. (lc) 

i <ij>crcr' 

Here i and j are site labels, the angle brackets in Eq. ( [la]) imply the restriction to neighboring 
sites, fii a = c irT Citj is the electron number operator on site i, and \i is the chemical potential. 
We will assume that the screened magnetic field inside the superconductor can be taken to 
be constant; for high T c materials this is a good approximation except for external fields 
extremely close to H c \. In a one-band lattice model the magnetic field appears in the hopping 
amplitudes: 

tij = te nc J n y ' (2) 

where t is the nearest neighbor hopping amplitude in zero field and V x A(r) = h(r). We 
report results below for two different models. The 's-wave' model has on-site attraction 
U < 0, and no nearest-neighbor interaction. For the 'd-wave' model, we set U > and 
V < 0. When the magnetic field is set to zero the mean-field BCS gap equations are readily 
solved for either model by using translational invariance on the lattice. For the 's-wave' case 
the pairing self-energy in the ordered state is proportional to the order parameter 



<3> s =< rLrt, > 



(3) 



while for the 'd-wave' model it is proportional to 



(4) 



where the sum is over the nearest-neighbors of site i, represented by the unit vectors 5 = 
±e x ,±e y . (In the homogeneous zero-field states $| and $f are independent of i.) For 
both models the numerical values of the interaction parameters have been chosen to give 
a zero-temperature coherence length (estimated from the pair wavefunction) ~ 4a, as in 
high T c materials. The results reported below were calculated for the case of a band filling 
factor < n >= 0.8. For the 's-wave' model, we set U = —3.5t and V = 0.0 while for the 
d-wave model we choose U = 2. It, and V = —2.1t. For layer separations and bandwidths 
appropriate for models of high-temperature superconductors, the penetration depths (at 
T = 0) for these models are several hundred times larger than the conherence lengths so 
that the models do indeed descripe strongly type-II superconductors. 

To study vortex lattice states, we introduce magnetic unit cells , each containing two 
superconducting flux quanta: $o = hc/2e. The size of a unit cell is N x a x N y a in general, 
where a is the lattice constant. We then define magnetic Bloch states labeled by a magnetic 
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wave vector k, a site index i within the magnetic unit cell, and a spin index a and denote 
the corresponding creation and annihilation operators by ct^ and cg icr . In mean-field theory 
these are related to the quasiparticle creation and annihilation operators by 



i t = E<(*)7L-EC(*)7& 4 , 



'feat 



(5a) 



-kii 



(5b) 



where u a (k) and v a (k) are determined by solving the BdG equations: 

H"{k) I /•!(/,•) t /•■>(/,•) //"(A') 

E a {k) 



H°(k) UF 1 (k) + VF 2 (k) 

UF*(k) +VF*(k) -H°*(k) J 
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v a (k) 



1 u a (k) ^ 



v a (k) 



(6) 



where H°, F\, and F 2 are N x N y x N x N y matrices. The off-diagonal blocks in this matrices are 
the pairing self-energies and these can be expressed in terms of quasiparticle amplititudes. 
The on-site interaction contribution is diagonal in site indices with 

mk)h =< c\^_ H >= -^tanh^ffi<(fc)<*(^), (7) 

while the nearest neighbor interaction contribution is 

(F 2 (k)) tJ =< cl/_ f . n > A j4+S = -W^nh^^-(uf(k)vf*(k) +u](k)vr(k))A j>i+s . 

(8) 

Here (3 = l/k B T and 

!Sj t i + s, if site i + 5 and site % are in the same cell; 
bj,i-N s s, if site i + 5 and site i are in different cells, 
where Ng is the number of sites along 5 direction in a cell. Eqs. (||), (0) and (H) constitute 
a set of self-consistent equations, whose solutions can be obtained numerically by iteration. 
Typical (|TJ| self-consistent results for the order parameter of the d-wave model in the 



vortex lattice state at T = are shown in Fig. [IJ. The size of the cell is 28a x 56a, corre- 
sponding to a field of H = $ /(28a) 2 ; if we associate a with the typical Cu-Cu distance in 
high T c materials this corresponds to a magnetic field ~ lOTesla. The value of the order 
parameter at the mid-point between neighboring vortices is ^f max {H) = 0.063, which may 
be compared with the value $ d = 0.065 obtained at T = for the Meissner state of the same 
model. The magnetic field suppresses superconductivity everywhere in the system; numeri- 
cal calculations at stronger fields are consistent |L5| with upper critical fields H c2 ~ lOOTesla 



for the model we study, as expected from the zero-temperature coherence length. We re- 
mark that the extended s-wave component of the order paramater, permitted by symmetry 
P JT2]JI3[] in the vortex lattice state and possibly of experimental relevance |T6| , |T7H in high T c 



materials, is smaller than the d-wave component by about two orders of magnitude for the 
model parameters we have chosen. 

We define the local density-of-states (LDOS) on site i by 



Ni(E) = -±^[\ut(k)\ 2 f'mk)-E) + \v?(k)\ 2 f'(E a (k) +E)l (9) 

k,a 

where f(E) = (exp(/3E) + and N c is the number of the magnetic cells in the system. 
Ni(E) is proportional to the differential tunneling conductance (§,[Tj| which is measured in 
an STM experiment. 

We show in Fig. Ni(E) at a site midway between two neighboring vortices for both 
s-wave and d-wave superconductors at zero field and at several different finite field strengths. 
In a field, the density of states at the Fermi level, iVj(O), is much larger for the d-wave case 
than for the s-wave case as predicted by Volovik. (We are unable to solve the BdG equations 
at weak enough fields to verify the expected H 1 ^ 2 behavior, although as seen in the inset our 
results are consistent with this prediction.) It is presumably this density-of-states away from 
the d-wave vortex cores which is responsible for the enhancement of the low-temperature 
specific heat of high T c superconductors in a field seen by Moler et al. [§,^(J. In the s-wave 
case, the zero-field gap remains quite well defined out to fairly strong magnetic field strengths 
although, in contrast with the semiclassical result, the density of states is not strictly zero 
at any energy. The size of the gap decreases with increasing field as expected. The sharp 
peaks in zero field DOS are of different origins. The peak closest to the Fermi energy is due 
to superconductivity while the second peak reflects the Van Hove singularity in the band 
structure. These peaks are smeared out in finite field. In Fig. |3] we show Ni(E) for a site at 
the center of a vortex core. In both s-wave and d-wave cases, we find large peaks near the 
Fermi energy, reflecting resonances which will evolve into quasiparticle bound states in the 
limit of isolated vortices. The positions of the peaks are distinctly different in two cases. 
In the d-wave model, the LDOS peak is not as strong and is clearly centered at the Fermi 



energy In the s-wave case two quasiparticle LDOS peaks are visible and the lowest 

energy of these is clearly located away from the Fermi energy. The scale of the separation 
between s-wave bound-state peaks is ~ 0.2t, in accord with expectations based on the size 
of the gap and the band width. For these short coherence length models, the separation 
between quasiparticle bound states is large enough to be resolved so that low-temperature 
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STM experiments would see a double-peak structure in high T c materials if they had s-wave 
symmetry, rather than the zero-bias peak observed || in conventional superconductors. 
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FIGURES 

FIG. 1. The amplitude of a d-wave order parameter &f (z axis) in a unit cell for the square 
lattice solution at T = 0. The size of the cell is 28a x 56a, corresponding to a field H = $o/(28a) 2 . 

FIG. 2. Quasiparticle local density-of-states profiles at T = away from vortex cores for 
different magnetic fields are plotted in (a) for the d-wave model and in (b) for the s-wave model. 
The largest cell size corresponds to the weakest magnetic field. The zero field limits for the same 
models are also plotted. All energies are in units of t and measured from the Fermi energy. The 
inset is the spatially averaged DOS N(0) vs H 1 / 2 . 

FIG. 3. Quasiparticle local density of states profiles at T = at the center of the vortex cores. 
Energies are in units of t and measured from the Fermi energy. 
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